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ABSTRACT 

The Wyoming Survey for Ha, or WySH, is a large-area, ground-based imaging 
survey for Ha-emitting galaxies at redshifts of z ~ 0.16, 0.24, 0.32, and 0.40. 
The survey spans up to four square degrees in a set of fields of low Galactic 
cirrus emission, using twin narrowband filters at each epoch for improved stellar 
continuum subtraction. Ha luminosity functions are presented for each Az ~ 
0.02 epoch based on a total of nearly 1200 galaxies. These data clearly show 
an evolution with lookback time in the volume-averaged cosmic star formation 
rate. Integrals of Schechter fits to the incompleteness- and extinction-corrected 
Ha luminosity functions indicate star formation rates per co-moving volume of 
0.010, 0.013, 0.020, 0.022 hm Mq yr'^ Mpc'^ at z ~ 0.16, 0.24, 0.32, and 0.40, 
respectively. Statistical and systematic measurement uncertainties combined are 
on the order of 25% while the effects of cosmic variance are at the 20% level. The 
bulk of this evolution is driven by changes in the characteristic luminosity of 
the Ha luminosity functions, with L^, for the earlier two epochs being a factor 
of two larger than at the latter two epochs; it is more difficult with this data 
set to decipher systematic evolutionary differences in the luminosity function 
amplitude and faint-end slope. Coupling these results with a comprehensive 
compilation of results from the literature on emission line surveys, the evolution 
in the cosmic star formation rate density over < 2; < 1.5 is measured to be 
Psfr(^)=Psfr(0)(1 + z)3-^±°-^ 
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Subject headings: galaxies: evolution — galaxies: luminosity function, mass func- 
tion 



1. Introduction 



The rate at which stars form on a cosmic scale is a fundamental property of the universe. 
Various efforts have shown the cosmic star formation rate per unit co-moving volume to 
significantly rise from the current epoch back to a redshift of 1, evolving at the rate of 
Psfr(-2) oc (1 + z)"^ , where 7 is estimated to lie between 1.5 and 4 (e.g., Lilly et al. 1996; 
Cowic ct al. 1999; Wilson et al. 2002; Gallego ct al. 2002; Trcssc ct al. 2002; Hippelein et al. 
2003; Hopkins 2004; Perez-Gonzalez et al. 2005; Le Floc'h et al. 2005; Hanish et al. 2006; 
Babbedge et al. 2006; Villar et al. 2008; Magnelh et al. 2009). Rehably constraining 7 will 
assist in interpreting models of galaxy evolution, since the cosmic star formation history 
is intimately tied to the build-up of heavy metals and stellar mass over time (Pei & Fall 
1995; Madau et al. 1996; Pei et al. 1999; Somerville et al. 2001; Cole et al. 2001; Panter 
et al. 2003). In fact, predictions of the accumulated stellar mass based on the integrated 
cosmic star formation history exceed the observed stellar mass at the present epoch by a 
factor of two (Figure 1 of Wilkins 2008). From an observational viewpoint, the evolutionary 
parameter 7 is ideally measured by uniformly sampling the star formation rate using a large 
sample of galaxies spanning multiple fields and multiple epochs between z = and ^ ~ 1 — 2. 
The Wyoming Survey for Ha, or WySH, surveys hundreds of galaxies at four intermediate 
redshifts, covering three separate fields that total some four square degrees (Dale et al. 
2008; hereafter D08). WySH utihzes Ha to probe the star formation in galaxies, an optical 
emission line that directly traces massive star formation, is technically simple to observe, 
and lies at a relatively long wavelength to help minimize the effects of extinction internal to 
star-forming galaxies. The survey has been carried out where redshifted Ha can be detected 
at wavelengths over which CCDs are sensitive: z ^ 0.16, 0.24, 0.32, and 0.40. This optical 
observing program is complemented by NewHa, a collaborative near-infrared narrowband 
imaging survey of Ha at redshifts oi z ~ 0.81 and 2.2 using the wide-field NEWFIRM 
camera (Probst et al. 2004) on the Kitt Peak National Observatory 4 m telescope (Lee et 
al. 2010; Ly et al. 2010). The NewHa survey includes observations through low airglow 
windows at 1.19 and 2.10 /im, sampling the Ha luminosity function down to levels similar 
to that explored in WySH. NewHa also has an extensive follow-up multi-slit spectroscopy 
campaign with the 6.5 m Baade telescope at Las Campanas Observatory, allowing us to 
verify Ha-emitting candidates, directly measure internal extinction via Balmer line ratios. 
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probe metal abundances, and to identify AGN (Momcheva et al. 2010). The combination of 
these two surveys will provide powerful constraints on the evolution in the Ha luminosity 
function and the cosmic star formation rate over < 2; < 2.2. This paper summarizes results 
from the WySH observing program, including Ha luminosity functions at four intermediate 
redshifts, a measure of the evolution in the volume-averaged cosmic star formation rate over 
these same epochs, and an estimate of the impact of cosmic variance. Section [2] reviews the 
survey parameters, § |3] presents results from all epochs, and § H] provides a summary. We 
assume Hq = 70 hjQ km s~^ Mpc^^, fl^a = 0.3, and fl\ = 0.7. 



2. The Survey 

WySH is a multi-epoch narrowband imaging survey for redshifted Ha carried out on the 
Wyoming Infrared Observatory 2.3 m telescope (WIRO). The survey spans multiple "blank" 
fields selected to minimize foreground contamination by zodiacal and Galactic dust, bright 
stars, and nearby galaxies. Observing multiple fields also aids in minimizing the impact of 
cosmic variance, a catchall term for the statistical fluctuations inherent to observations of 
different regions at cosmological distances. By design, two of these fields, ELAIS-Nl and 
the Lockman Hole, overlap with the target areas of deep infrared and ultraviolet surveys. 
The SWIRE Spitzer Legacy project provides maps of these regions using several bandpasses 
between 3 and 160 /xm (Lonsdale et al. 2003) along with a wealth of ancillary data at other 
wavelengths. The GALEX Deep Imaging Survey (Martin et al. 2005) includes the ELAIS-Nl 
and Lockman Hole fields, providing two channels of ultraviolet data from integrations that 
are 300 times longer than the standard integration of that mission's all-sky survey. The third 
WySH field, also known as "WySH 1", conveniently fills the right ascension gap between the 
Lockman Hole and the ELAIS-Nl field to enable year-round observations from WIRO (DOS 
Table 1). While the main goal of this survey is to provide a statistically robust measure of 
the evolving star formation history of the Universe at intermediate redshifts, the combination 
of the surveys described above enables interesting parallel science. For example, combining 
Ha, ultraviolet, and infrared data allow us to study the evolution over time of the average 
attenuation by interstellar dust within star-forming galaxies (Moore et al. 2010). 



2.1. Observations and Data Processing 

WIRO is equipped with a prime focus CCD camera with 0'.'523 pixels and a 17^9 field 
of view (Pierce & Nations 2002). Redshifted Ha emission is observed in four separate 
cosmological epochs, z ~ 0.16,0.24,0.32, and 0.40. This survey differs from traditional 
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narrowband imaging surveys by using nearly identical and wavelength- adjacent narrowband 
filter pairs (~ 60A FWHM). If a source's redsliift places the Ha emission to appear in 
Filter A, then the emission detected from Filter B is used to subtract off the stellar continuum 
as detected by Filter A. Likewise, the reverse applies if Ha is found in Filter B. The basic 
unit of observation is a 300 s frame, though several are taken for each filter at each location. 
After the pre-processing of individual frames (e.g., bias subtraction, flat- fielding, fringing 
correction; see § 2.3 of DOS), the multiple 300 s frames for a given field are ahgned and 
stacked to create images with longer effective integrations. The effective integration increases 
with increasing redshift with the aim of achieving a similar sensitivity to line luminosity at 
each epoch: the z Ri[0.16, 0.24, 0.32, 0.40] stacks have effective integrations of [1200, 3600, 
6000, 9600] s. Since two filters are utilized at each epoch, there are in fact two stacks of 
these integration depths for each pointing at each epoch. The total areal coverage for the 
survey is [4.19, 4.03, 4.13, 1.11] square degrees at z fti[0.16, 0.24, 0.32, 0.40]. 

The details of how sources are identified and their photometry extracted with SExtractor 
(Bertin & Arnouts 1996) is described in Dale et al. (2008). Salient details include a 5% flux 
calibration uncertainty, a 20% correction for [N II] emission within the filter bandpasses, 
and the application of a 3o" cut based on the Ha+[N II] signal-to- noise. The 3o" continuum 
sensitivity to luminosity is ~ 10^° ergs s~^ (and ~ 5 ■ 10~^^ ergs s~^ cm~^ in terms of fiux), 
with the survey at the later epochs being somewhat more sensitive and the surveys at the 
earlier epochs somewhat less sensitive than this value, in part due to the falling quantum 
efficiency of the CCD with wavelength (Figure 5 of DOS). The sensitivity to line detections 
is an additional factor of ~ -\/2 poorer since line strengths are derived from a subtraction of 
two narrowband luminosities. 



2.2. Removal of Contamination by Redshift Interlopers 

The inclusion of non-Ha line emission from galaxies outside of the target epochs will bias 
the Ha luminosity function and thus the inferred star formation rate density at each epoch; 
other prominent optical emission lines such as [O IIl]AA4959/5007, [O II]A3727, and H^A4S61 
could be rcdshifted into the filter bandpasses (e.g., Fujita et al. 2003). Photometric redshifts 
are used to cull the contaminators. Multiple photometric redshift catalogs are available for 
the fields pursued in this work: Rowan-Robinson et al. (2007) from SWIRE, Csabai et al. 
(2003) from SDSS, and from our own suite of deep observations from WIRO covering UBVRI 
(15-75 min integrations) and narrowband imaging at eight wavelengths between 7597 and 
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9233 nm. WIRO-based photometric redshifts are estimated using the Le Phare codqll and 
a suite of spiral, elhptical, irregular, and starburst templates (the CWW_KINNEY templates), 
allowing the V band extinction to vary from to 3 magnitudes. For the sake of uniformity 
SWIRE photometric redshifts are used where available, with SDSS and our own photometric 
redshifts employed to supplement the SWIRE coverage. 



3. Results 

3.1. Ha Luminosity Functions 

The Ha luminosity function for each epoch is displayed in Figure [H The amplitude of 
the i^^ bin of each luminosity function and its uncertainty are derived from 

where V{zj,Lj) is the co- moving volume for the j^^ galaxy in the summation. The lumi- 
nosities are binned according to | logLj — logLj | < ^AlogL, with a luminosity bin width 
AlogL spanning 0.4 dex. Several corrections described in § 3 of DOS are incorporated into 
the luminosity functions displayed in Figure [TJ The luminosity functions are corrected for 
sample incompleteness i^)inc), filter transmission characteristics, and Ha line extinction 
(^gT(i)^ following the Ha luminosity- dependent prescription of Hopkins et al. (2001; though 
see also Garn et al. 2010). A bin- by-bin tabulation of the incompleteness and extinction 
corrections are provided in Table [TJ Open circles in Figure [1] indicate the data corrected for 
all issues described above except incompleteness; the filled circles also include corrections 
for incompleteness. Error bars in Figure [T] refiect the uncertainty in the luminosity function 
amplitude according to Equation [H summed in quadrature with the uncertainties in the 
incompleteness corrections. Also included in Figure [1] are Ha-based luminosity functions 
from the literature, largely consistent with our results. 

Going down to and including these limiting luminosities, a Schechter profile (Schechter 
1976) is fit to the incompleteness-corrected luminosity functions: 

$(z,logL)rflogL = (t){z,L)dL = 0,(z) -— e-^/^*(^)-^, (2) 



where a{z) conveys the shape of the function, L*(z) sets the luminosity scale, and </)*(-2) 
represents the overall normalization. The parameters for the functional fits displayed in 



^http://www. oamp.fr/people/arnouts/LE_P HARE. html 
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Figure [T] are listed in the first four rows of Table [2l Since the faint end slopes are difficult 
to gauge from the available luminosity bins for the two earliest epochs, for these solutions 
we have fixed a{z ~ 0.32) and a{z ^ 0.40) to be the same as the average of that found for 
the later two epochs (—1.38). Also included in Table [2] are solutions allowing all parameters 
to vary. In all cases Monte Carlo simulations are employed to estimate the average fit 
parameters and their uncertainties (see Figure |2]). Each fit is repeated 10,000 times after 
using the measured luminosity uncertainties to randomly add a (Gaussian deviate) offset to 
each Ha luminosity. The standard deviations in the simulations are used to represent the 
luminosity function parameter uncertainties. 

Extremely deep surveys are required to adequately measure the faint end slopes of the lu- 
minosity function (e.g., Ly et al. 2007). As can be seen from Figure [H the lowest luminosities 
for which ^{z, log L) are reliably obtained increase with redshift: log Liijn(Ha/ergs s^^) =[40.2, 
40.6, 41.0, 41.4] for z ^ [0.16, 0.24, 0.32, 0.40]. While the fitted values for a{z) are somewhat 
steeper for z ^ 0.16 and 0.24 than at the two earlier epochs (see Rows 5-8 in Table |2]), this 
result is very tenuous since there are too few low luminosity galaxies detected at z ^ 0.32 
and 0.40 to adequately constrain the faint end slopes. 

In terms of the other luminosity function parameters, L*(z) and </)*(-2), Westra et al. 
(2010) show that intermediate redshift surveys like WySH that span at least 3 square degrees 
cover large enough volumes to effectively constrain the bright end of the luminosity function. 
For our standard solutions involving a fixed slope of a = —1.38 and luminosity-dependent 
internal extinction, the changes in the luminosity function with redshift are mainly driven 
by changes in the characteristic luminosity L^,{z), with the average L^,(z) being a factor of 
two larger for the two earlier epochs than the average measured at 2; ~ 0.16 and 0.24. It is 
more difficult to discern any systematic evolutionary changes in the source number density 
as parameterized by the luminosity function amplitude 0*(2;). 



The volume-averaged cosmic star formation rate can be computed by integrating under 
the fitted Schechter function and multiplying by the Kennicutt (1998) star formation rate 
calibration: 



3.2. The Evolution in the Cosmic Star Formation Density 



PsFRiz)ihjo Mq yr-1 Mpc^^^) = 7.9 ■ 10^^^ C{z){hjo ergs Mpc-=^) 



(3) 



where an analytical expression for the luminosity density is 




(4) 
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Table |2] provides the integrated cosmic star formation rate densities. There is a clear, sys- 
tematic increase in Psfr(^) with redshift, with an overall difference of a factor of ~ 2 between 
z ^ 0.16 and z ^ 0.40. This result is robust to the effects of cosmic variance, fluctuations 
due to the characteristics of the particular volume(s) being probed along a survey's line(s)- 
of-sight (e.g., clusters, voids, etc.): the values for Psfr(^) extracted from individual fields 
(e.g., Lockman Hole and ELAIS-Nl) differ by only ~ 20%. Placed in a larger context, our 
values for Psfr{z) fit well within the envelope determined by previous emission line surveys 
(see Figure Eland the citations in the caption). If the evolution over < 2; < 1.5 is cast in 
power-law form, Psfr(-2) oc (1 + , the exponent is best fit by 7 = 3.4 ± 0.4. This slope is 
consistent with the typical values (3 < 7 < 4) found in the literature for this redshift range 
using star formation rates based on the infrared continuum or dust-corrected optical emis- 
sion lines (e.g., Gallego et al. 2002; Tresse et al. 2002; Hippelein et al. 2003; Perez-Gonzalez 
et al. 2005; Le Floc'h et al. 2005; Hanish et al. 2005; Babbedge et al. 2006; Villar et al. 
2008; Magnelli et al. 2009). Surveys based on the ultraviolet continuum, where attenuation 
by dust and contributions from older stellar populations can play important roles, tend to 
recover shallower slopes (e.g., Cowie et al. 1999; Wilson et al. 2002; Schiminovich et al. 
2005; Prescott et al. 2009). Cosmic star formation history surveys using optical emission 
line data at higher redshifts are understandably sparse, and Figure |3] provides just a hint of 
the leveling off of psfr{z) beyond a redshift of 2; ~ 1 — 1.5 that is seen in other compilations 
that profit from a broader wavelength baseline (e.g., Hopkins 2004). 



4. Discussion and Summary 



The primary aim of the Wyoming Survey for Ha, or WySH, is to accurately quantify the 
Ha luminosity function via narrowband imaging spanning < 4 square degrees and multiple 
cosmic epochs. Important features of the survey include the use of narrowband filter pairs 
at each epoch for improved stellar continuum subtraction and the spatial overlap with deep 
ultraviolet and infrared surveys that enable interesting follow-up studies. Buttressed by a 
total of nearly 1200 Ha detections, we find a modest evolution in the Ha luminosity function 
over z ~ 0.16, 0.24, 0.32, and 0.40. The values of the volume-averaged cosmic star formation 
rate, found by integrating under the luminosity functions, change by a factor of two over 
this moderate stretch in redshift. Our results indicate that this evolution is largely driven 
by changes in the characteristic luminosity L*, which also shows an evolution by a factor of 
two over these epochs. That the evolution in the cosmic star formation rate density over 
these intermediate redshifts is mainly influenced by systematic changes in the characteristic 
luminosity is consistent with the findings of Le Floc'h et al. (2005), Perez- Gonzalez et al. 
(2005), MagneUi et al. (2009), and Westra et al. (2010), though Ly et al. (2007) find the 
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evolution to be more driven by changes in the source number density. Placing our results in 

the larger context of the slew of recent emission line surveys for Psfr over < 2; < 1.5, the 
evolution in the cosmic star formation rate density is estimated to be Psfr{z) oc (l + z)^'*^^"'"^. 
Results from a complementary near-infrared narrowband imaging survey of H«-emitters will 
extend this work to redshifts of 2; ~ 0.81 and 2.2. Finally, the large volume covered by this 
optical survey enables a measure of the impact of cosmic variance. By separately analyzing 
the different fields in this survey, in particular the Lockman Hole and ELAIS-Nl, we find a 
variation in psfr at the 20% level for any given redshift. 

This research is funded through the NSF CAREER program (AST0348990). Special 
thanks go to Michael Pierce for developing WIROPrime and to the WIRO staff for keeping 
the telescope operational. IRAF, the Image Reduction and Analysis Facility, has been de- 
veloped by the National Optical Astronomy Observatories and the Space Telescope Science 
Institute. 
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Table 1. Ha Extinction and Luminosity Function Incompleteness 



logL(Ha) 






l^i^t ^)inc 


^)inc 


^)inc 


(ergs s~i) 




7597/7661A 


8132/8199A 


8615/8685A 


9155/9233A 


40.2 


1.18 


0.15±0.02 








40.6 


1.49 


0.65±0.05 


0.27±0.04 






41.0 


1.84 


0.89±0.05 


0.66±0.07 


0.15±0.05 




41.4 


2.24 


0.91±0.04 


0.89±0.05 


0.69±0.05 


0.21±0.05 


41.8 


2.68 


0.93±0.04 


0.92±0.04 


0.91±0.04 


0.79±0.04 


42.2 


3.17 


0.96±0.03 


0.96±0.03 


0.96±0.03 


0.95±0.03 


42.6 


3.70 






0.98±0.03 


0.96±0.03 



Table 2. Luminosity Function Results 



Redshift 


a 


logL* 


log0. 


PSFR 






{hjo ergs s-i) 


(^gMpc-^dex-i) 


/i7oM0yr-lMpc-3 


0.16 


-1.38^ 


42.0±0.1 


-3.08±0.02 


O.OlOiO.OOl 


0.24 


-1.38'' 


41.8±0.1 


-2.70±0.02 


0.013±0.001 


0.32 


-1.38" 


42.2±0.1 


-2.90±0.01 


0.021±0.001 


0.40 


-1.38" 


42.3±0.1 


-2.97±0.04 


0.024±0.002 


0.16 


-1.36±0.06 


42.0±0.2 


-3.05±0.11 


O.OlOitO.OOl 


0.24 


-1.41±0.05 


41.8±0.1 


-2.74±0.07 


0.013±0.001 


0..'?2 


-1.26±0.0.5 


42.1±0.1 


-2.77±0.0.'5 


0.020±0.001 


0.10 


-1.11=0.21) 


12.2±().l 


-2,7<)=0,l(i 


0.022±0.002 


''Fixed; 


not fitted. 
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I — 1 Star Formation Rate (Mq yr"^) 

^ 0.1 1 10 0.1 1 10 0.1 1 10 0.1 1 10 




41 42 41 42 41 42 41 42 



log[L(Ha)/ergs s"^] corrected for extinction 



Fig. 1. — The luminosity functions at z ~ 0.16, 0.24, 0.32, and 0.40, based on 
214+424+438+91 Ha-emitting galaxies. The data without incompleteness corrections are 
displayed as open circles, while those corrected for incompleteness are shown as filled cir- 
cles. Error bars reflect the uncertainty in the luminosity function amplitude according to 
Equation [H summed in quadrature with the uncertainties in the incompleteness corrections. 
The thick solid lines show the Schechter fits for the parameters presented in the first four 
rows of Table O Literature fits are also provided at 2; ~ 0.16 (Sullivan et al. 2000; Westra 
et al. 2010), z ~ 0.24 (Tresse & Maddox 1998; Fujita et al. 2003; Ly et al. 2007; Shioya et 
al. 2008), and z ~ 0.40 (Ly et al. 2007). 
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Fig. 2. — The 2cr (95.4%) confidence distributions of tlie luminosity function parameters a 
and for 10,000 Monte Carlo simulations of the data, for the case where all parameters are 
allowed to vary (Rows 5-8 of Table [2]). The redshifts are indicated within the contours. 
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Fig. 3. — The cosmic star formation history, using results from hne emission surveys only. 
Data points from the WySH survey presented here are encircled and have a slope of 4.5±0.7. 
Literature data come from Gallego et al. (1995), Hogg et al. (1998), Tresse & Maddox (1998), 
Yan et al. (1999) assuming Ay = I mag, Sullivan et al. (2000), Gallego et al. (2002), Tresse 
et al. (2002), Fujita et al. (2003), Perez-Gonzalez et al. (2003), Hippelein et al. (2003), 
Nakamura et al. (2004), Brinchmann et al. (2004), Hanish et al. (2006), Takahashi et al. 
(2007), Ly et al. (2007), Geach et al. (2008), Westra & Jones (2008), Morioka et al. (2008), 
Shioya et al. (2008), Villar at al. (2008), Dale et al. (2008), Shim et al. (2009), Sobral et al. 
(2009), Ly et al. (2010), Hayes et al. (2010), and Westra et al. (2010). 



